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We propose quasiperiodic iieterostructures associated with the tessellations of the unit disk by regular hyper- 
bolic triangles. We present explicit construction rules and explore some of the properties exhibited by these 
geometric-based systems. 
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Quasiperiodic (QP) systems have been receiving a lot of at- 
tention over the last years \V\. The interest was originally mo- 
tivated by the theoretical predictions that they should manifest 
peculiar electron and phonon critical states 121 |3|l, associated 
with highly fragmented fractal energy spectra ij, S |6[1 . On 
the other hand, the practical fabrication of Fibonacci ||7|] and 
Thue-Morse [8] superlattices has triggered a number of exper- 
imental achievements that have provided new insights into the 
capabilities of QP structures 191]. In particular, possible optical 
applications have deserved major attention and some intrigu- 
ing properties have been demonstrated il0lllUll2lll3ll . Un- 
derlying all these theoretical and experimental efforts a cru- 
cial fundamental question remains concerning whether QP de- 
vices would achieve better performance than usual periodic 
ones for some specific applications il4ll . 

The QP systems considered thus far rely for their explicit 
construction in substitutional rules among the elements of a 
basic alphabet. In the common case of a two-letter alpha- 
bet {A.B}, the algorithm takes the form A h-> <ja{A,B), 
B I— > (Jb{A, B), where a a and as can be any string of the 
letters. The sequences generated after n applications of the 
algorithm are of significance in fields as diverse as cryptogra- 
phy, time-series analysis, and cellular automata fl5]. In ad- 
dition, they have interesting algebraic properties, which are 
usually characterized by the nature of their Fourier or multi- 
fractal spectra jldl . 

We wish to approach the problem from an alternative ge- 
ometrical perspective. To this end, we first observe that in 
many problems of physical interest OlTIl the letters of the al- 
phabet can be identified with one-dimensional linear lossless 
systems (i.e., with two input and two output channels). Un- 
der these general conditions, it turns out that the associated 
transfer matrix belongs to the group SU(1,1), which is also 
the basic symmetry group of the hyperbolic geometry 1 18]. In 
consequence, the unit disk a ppe ars as the natural arena to dis- 
cuss their performance ll l9[ \2(JL 12111 . Since in the Euclidean 
plane, QP behavior is intimately linked with tessellations, one 
is unfailingly led to consider the role of hyperbolic tessella- 
tions in the unit disk, much in the spirit of Escher's master- 
piece woodcut Circle Limit III ll22ll . The answer we propose 
is promising: the tessellations by different regular polygons 
provide new sequences with properties that may open avenues 



of research in this field. 

A QP system can thus be seen as a word generated by 
stacking different letters of the basic alphabet. To be spe- 
cific, we focus our attention on the optical response. Let us 
consider one of these letters (which in practice is made of 
several plane-parallel layers), which we assume to be sand- 
wiched between two semi-infinite identical ambient (a) and 
substrate (s). We suppose monochromatic plane waves in- 
cident, in general, from both the ambient and the substrate. 
As a result of multiple reflections in all the interfaces, the to- 
tal electric field can be decomposed in terms of forward- and 
backward-traveling plane waves, denoted by _E(+) and £'(^\ 
respectively. If we take these components as a vector 
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then the amplitudes at both the ambient and the substrate sides 
are related by the transfer matrix M 
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ME, 



It can be shown that M is of the form 
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where the complex numbers R and T are, respectively, the 
overall reflection and transmission coefficients for a wave in- 
cident from the ambient. The condition dot M = +1 is equiv- 
alent to |i?p + |Tp ~ 1, and then the set of transfer matrices 
reduces to the group SU(1,1). Obviously, the matrix of a word 
obtained by putting together letters of the alphabet is the prod- 
uct of the matrices representing each one of them, taken in the 
appropriate order. 

In many instances we are interested in the transformation 
properties of field quotients rather than the fields themselves. 
Therefore, it seems natural to consider the complex numbers 



E(- 



Ei+) 
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for both ambient and substrate. The action of the matrix given 
in Eq. ^ can be then seen as a function Za — f{zs) that can 
be appropriately called the transfer function. From a geomet- 
rical viewpoint, this function defines a transformation of the 



complex plane C, mapping the point Zs into the point Za ac- 
cording to 
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where a — l/T and fi — R* /T*. When no Hght is incident 
from the substrate, Zg — and then Za = R- Equation (|5]) is 
a bilinear (or Mobius) transformation. One can check that the 
unit disk, the external region and the unit circle remain invari- 
ant under (|5]l. This unit disk is then a model for hyperbolic 
geometry in which a line is represented as an arc of a circle 
that meets the boundary of the disk at right angles to it (and 
diameters are also permitted). In this model, we have three 
different kinds of lines: intersecting, parallel (they intersect 
at infinity, which is precisely the boundary of the disk) and 
ultraparallel (they are neither intersecting nor parallel). 

To classify the possible actions it proves convenient to work 
out the fixed points of the transfer function; that is, the field 
configurations such that Za — Zs = Zf in Eq. Q, whose 
solutions are 



Zf 



2/3 I 



2i Im(a) ± A/[Tr(M)]2 - 4| 
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When [Tr(M)]^ < 4 the action is elliptic and it has only one 
fixed point inside the unit disk. Since in the Euclidean geom- 
etry a rotation is characterized for having only one invariant 
point, this action can be appropriately called a hyperbolic ro- 
tation. 

When [Tr(M)]^ > 4 the action is hyperbolic and it has 
two fixed points, both on the boundary of the unit disk. The 
geodesic line joining these two fixed points remains invariant 
and thus, by analogy with the Euclidean case, this action is 
called a hyperbolic translation. 

Finally, when [Tr(M)]^ — 4 the action is parabolic and it 
has only one (double) fixed point on the boundary of the unit 
disk. 

As it is well known, SU(1, 1) is isomorphic to the group of 
real unimodular matrices SL(2, M), which allows us to trans- 
late the geometrical structure defined in the unit disk to the 
complex upper semiplane, recovering in this way an alterna- 
tive model of the hyperbolic geometry that is useful in some 
applications. 

The notion of periodicity is intimately connected with tes- 
sellations, i.e., tilings by identical replicas of a unit cell (or 
fundamental domain) that fill the plane with no overlaps and 
no gaps. Of special interest is the case when the primitive cell 
is a regular polygon with a finite area |23|]. In the Euclidean 
plane, the associated regular tessellation is generically noted 
{p, q], where p is the number of polygon edges and q is the 
number of polygons that meet at a vertex. We recall that ge- 
ometrical constraints limit the possible regular tilings {p, g} 
to those verifying {p — 2){q — 2) = 4. This includes the 
classical tilings {4, 4} (tiling by squares) and {6, 3} (tiling 
by hexagons), plus a third one, the tiling {3, 6} by triangles 
(which is dual to the {6, 3}). 
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FIG. 1: (color online). Arealisticimplementationof the generators A 
and B (together with A^^) of the tessellations by hyperbolic triangles 
in (O. The external layers (in blue) are made of cryolite (NasAlFs), 
while the central medium (in brown) is zinc selenide (ZnSe). The 
wavelength in vacuum is A = 610 nm and normal incidence (from 
left to right) has been assumed. The corresponding thicknesses are 
expressed in nanometers. 



On the contrary, in the hyperbolic disk regular tilings exist 
provided [p — 2)(g — 2) > 4, which now leads to an infinite 
number of possibilities. An essential ingredient is the way 
to obtain fundamental polygons. These polygons are directly 
connected to the discrete subgroups of isometrics (or congru- 
ent mappings). Such groups are called Fuchsian groups ll24ll 
and play for the hyperbolic geometry a role similar to that of 
crystallographic groups for the Euclidean geometry 12511 . 

A tessellation of the hyperbolic plane by regular poly- 
gons has a symmetry group that is generated by reflections in 
geodesies, which are inversions across circles in the unit disk. 
These geodesies correspond to edges or axes of symmetry of 
the polygons. Therefore, to construct a tessellation of the unit 
disk one just has to built one tile and to duplicate it by using 
reflections in the edges. 

In this Letter, we consider only the simple example of a tes- 
sellation by triangles with vertices in the unit circle, although 
the treatment can be extended to other polygons. The key 
idea is to consider the Fuchsian group generated by an ellip- 
tic transformation whose fixed point is the middle of an edge 
of the triangle and a parabolic one with its fixed point in the 
opposite vertex. Proceeding in this way we get 



l + i/Vi 1/V3 
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The fixed point of A is — i; while for B the fixed point in the 
disk is i(2 — \/3). In Fig. 1 we show a possible way in which 
these matrices can be implemented in terms of two commonly 
employed materials in optics. Note that, in physical terms, 
the inverses must be constructed as independent systems, al- 
though in our case only A^^ must be considered, since the 
action of B^^ coincides with that of B. 

In Fig. 2 we have shown the tessellation obtained by trans- 
forming the fundamental triangle with the Fuchsian group 
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FIG. 2: (color online). Tiling of the unit disk with the matrices ([TJ. 
The marked points are the barycenters of the triangles in the tessel- 
lation and all of them are the transformed of the origin by a matrix 
that have as reflection coefficient the complex number that links the 
origin with the center of the triangle. 



generated by the powers of {A, B} (and the inverses). This tri- 
angle is equilateral with vertices at the points — i, exp(J7r/6) 
and cxp(i57r/6) (which are the fixed points of A, AB and BA, 
respectively). Moreover, all the other triangles are equal, with 
an area tt. In the figure we have plotted also the barycenters of 
each triangle together with the resulting tree (that is called the 
dual graph of the tessellation), which turns out to be a Farey 
tree 12611 . In fact, each line connecting two of these barycen- 
ters represent the action of a word (with alphabet {A, B} and 
the inverses). 

To give an explicit construction rule for the possible words, 
we proceed as follows. First, we arbitrarily assign the num- 
ber to the upper side of the fundamental triangle, while the 
other two sides are clockwise numbered as 1 and 2. It is easy 
to convince oneself that this assignment fixes once for all the 
numbering for the sides of the other triangles in the tessella- 
tion. However, these triangles can be distinguished by their 
orientation (as seen from the coiTesponding barycenter): the 
clockwise oriented are filled in red, while the counterclock- 



TABLE I: Explicit rules to obtain the barycenter Zn+i from the z„. 
We have indicated the corresponding transformations, which depend 
on the color jumps and the sides crossed by going from z„ to z„+i. 
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FIG. 3: (color online). Normalized structure factor for the word 
(made of 60 letters) connecting the origin with the point za in the 
zig-zag path shown in Fig. 2. 



wise are filled in yellow. In short, we have determined a fun- 
damental coloring of the tessellation Ii27i] . 

To obtain one barycenter Zn+i from the previous one Zn, 
one looks first at the corresponding color jump. Next, the ma- 
trix that take z„ into Zn+i depends on the numbering of the 
side (0, 1, or 2) one must cross, and appear in the appropriate 
column "Trans" in Table |T] The next generation is obtained 
much in the same way, except for the fact that A„ and B„ must 
be replaced by A„+i and B„+i, respectively, as indicated in 
the Table. In obtaining recursively any word, the origin is de- 
noted as zo and the matrices Aq and Bq coincide with A and 
B. 

With this rule, one can construct any word proceeding step 
by step. For example, the word that transform zq into z^ in the 
zig-zag path sketched in Fig. 2 results 



(8) 



Obviously, the total word is obtained by composing these par- 
tial words. 

In fact, one can show that given an arbitrary sequence of 
nonzero integers {ki, fc2, . . . , A:,}, the word represented by 
the transfer matrix 
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where {si, S2} C {0, 1}, transforms the origin in a barycenter 
of the tessellation. 

Given the geometric regularity of the construction sketched 
in this Letter, the sequences obtained must play a key role in 
the theory and practice of QP systems. Of course, to put for- 
ward the relevant physical features of these sequences, there 
are a number of quantities one can look at. Perhaps one of 
the most appropriate ones to assess the performance of these 
systems is the structure factor [15]. For a word M (composed 
of L letters) we define a numerical sequence fn by assigning 
1, e2iri/3^ ^jj(j g-27ri/3 f^ ^^ letters B, A, and A~\ respec- 
tively. Next, we calculate the discrete Fourier transform of the 
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where k — 0,1, ..L — 1. The structure factor (or power spec- 
trum) is just |Ffcp. In Fig. 3 we have plotted this structure 
factor in terms of k for the word connecting the origin with 
the point zg in the zig-zag path of Fig. 2. The peaks reveal a 
rich behavior: a full analysis of these questions is outside the 
scope of this Letter and will be presented elsewhere. 

As a final and rather technical remark, we note that the quo- 
tient of the hyperbolic disk by the Fuchsian group generated 
by A and B is a 2-orbifold of genus 0, with a conical point of 
order two and a cusp. Each word as given in Eq. (|9]l repre- 
sents a hyperbolic transformation of the disk, and the axis of 
the transformation is projected onto a closed geodesic of such 
an orbifold. This provides an orbifold interpretation of our QP 
sequences. 

In summary, we expect to have presented new schemes to 
generate QP sequences based on hyperbolic tessellations of 
the unit disk. Apart from the intrinsic beauty of the formalism, 
our preliminary results seem to be quite encouraging for future 
applications of these systems. 

The authors wish to express their warmest gratitude to E. 
Macia and J. M. Montesinos for their help and interest in the 
present work. 
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